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We discuss a five-dimensional Minimal Supersymmetric Standard Model compactified on a S /Z% 
orbifold, looking at, in particular, the one-loop evolution equations of the Yukawa couplings for the 
quark sector and various flavor observables. Different possibilities for the matter fields are discussed, 
that is, where they are in the bulk or localised to the brane. The two possibilities give rise to quite 
1 different behaviours. By studying the implications of the evolution with the renormalisation group 

of the Yukawa couplings and of the flavor observables we find that, for a theory that is valid up to 
the unification scale, the case where fields are localised to the brane, with a large tan ft, would be 
£SJ 1 more easily distinguishable from other scenarios. 

I. INTRODUCTION 

The masses of the quarks and charged leptons are determined in the Standard Model (SM) via Yukawa couplings 
to the Higgs boson. The origin of their structure (masses and mixing angles) has no explanation within the SM and 
Q-i presents as one of the major challenges for physics beyond the SM. Among these models those with extra spatial 
^ ■ dimensions offer many possibilities for model building and TeV-scale physics scenarios which can now be explored or 
constrained at the Large Hadron Collider (LHC). Extra-dimensional models allow, for example, a way to generate 
. electroweak symmetry breaking or supersymmetry breaking through the choice of appropriate boundary conditions 
(for a review see Ref.[l|). In addition, for the case of flat extra-dimensions, the presence of towers of excited Kaluza- 
. Klein (KK) states induces a power-law enhancement of the gauge couplings, leading to possible low-scale unification 
0, This effect can be applied to other couplings, such as Yukawa couplings, giving an original way to generate 
£NJ [ mass hierarchies [U, Q . The study of the renormalisation group equations (RGEs) provides a way by which partial 
■ explorations of the physics implications at a high energy scale is possible, as the theories at asymptotic energies may 
reveal new symmetries or other interesting properties which may lead to deeper insights into the physical content. As 
such, in order to understand these as yet unexplained hierarchies of mixing angles and fermion masses, a great deal 
of work has gone into analysing the RGEs |5l-[l2|. Therefore the behaviour of the well known quark sector's flavor 
mixings in the charged current, as described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix, shall be a focus of 
this paper. 

By using the RGEs, we will study the asymptotic properties of the Lagrangian parameters like the Yukawa coupling 
constants and mixing angles, and explore the possibility of a model in which the CKM matrix might have a simple, 
special form at asymptotic energies. This quest will, however, be limited by the fact that extra-dimensional theories 
are only effective ones, limited by a cut-off in their physical description of fundamental phenomena. Therefore 
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the following study can only be used as an indication of the behaviour of couplings and mixing parameters at an 
intermediate scale between the electroweak scale (at which these parameters are measured) and the higher scale at 
which the effective theory ceases to be valid. However, in the range of the LHC energies (of the order of a few TeV), 
and beyond, one can indeed test if the departure from the usual behaviour of the coupling evolution can be seen in 
precision flavor measurements. In order to discuss the implications of an effective five dimensional supersymmetric 
theory, we first recall the usual formalism in the SM, where the CKM matrix has four observable parameters including 
three mixing angles and one phase. In the standard parametrisation it has the form: 

V us V ub \ I C12C13 S12C13 

V C KM = I V cd V cs V cb = -S12C23 - Ci 2 S 2 3Sl3e l<5 C12C23 - S 12 S 2 3S 13 e lS S 23 C 13 j (1) 

Vtd V ts V tb J \ s 12 s 23 - c 12 c 23 Si 3 e' lS -C12S23 - si 2 c 23 s 13 e lS 
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where S12 = sin #12, cyi = cos #12 etc. are the sines and cosines of the three mixing angles #12, #23 and #13, and 6 is 
the CP violating phase. 

The aim of this paper is to therefore study these effects explicitly in the case of one extra-dimension within a minimal 
supersymmetric model, which will be referred to as 5D MSSM. Note that as pointed out in Ref.[l3| a supersymmetric 
and extra-dimensional extension to the SM can be considered simultaneously, with the advantages that though higher 
dimensional field theories are non-renormalizable (requiring the existence of a UV completion), superstring theories 
(which contain supersymmetry) provide the main hope in this context. In fact, such a point of view has already 
been explored through AdS/CFT correspondences in accounting for both "little" and "big" hierarchies [14j . Also 
embedding the SM in an extra-dimensional space-time does not stabilise instabilities due to quantum corrections in 
the scalar potential, whilst supersymmetry may ameliorate this. Finally, most 4-dimensional supersymmetric models 
lack a simple mechanism for supersymmetry breaking, something the extra-dimensions may offer. To be precise, we 
will discuss a five dimensional Af — 1 supersymmetric model compactified on the S 1 /Z2 orbifold as a simple testing 
ground for the effects of the extra-dimension on the quark Yukawa couplings and the CKM matrix observables. We will 
not discuss in the present paper the RGE evolution of the sparticle mixing parameters sector. This interesting question 
will be explored elsewhere as it requires more detailed assumptions on the sparticle spectrum and corresponding soft 
supersymmetry breaking (SSB) parameters, where at least some SSB parameter matrices can be non-diagonal in the 
basis in which quarks are diagonal (this implies flavor- violating decays, giving rise to the well known flavor problem 
in the MSSM and its extensions; even without considering the RGE evolution of the SSB parameters). 

The paper is organised as follows: In section [TT] we shall develop our model and calculate the various beta functions 
required to determine the evolution equations for the Yukawa couplings. We shall look at two cases, where the matter 
fields are localised on the brane, and where they are free to propagate in the bulk. In section IIIII we shall then 
determine the evolutions of fermion Yukawa couplings, and the CKM matrix which follows from them. In section [IV] 
we then numerically analyse the evolution behaviours of the physical observables, including the Yukawa couplings, 
quark flavor mixings, and the CP violation Jarlskog parameter. Our results and conclusions are presented in section 



II. THE 5D MSSM AND BETA FUNCTIONS 



In this Universal Extra-dimensional (UED) formalism the compactification S1/Z2 implies that the 5D Lorentz 
invariance is broken to the usual 4D one, however, a remnant of momentum conservation along the fifth coordinate 
implies that KK number is conserved at tree level and that KK parity is conserved at loop level. The KK-parity 
invariance has two main implications in phenomenology: the contributions of the KK modes to electroweak precision 
observables arise only through loops, and the exact KK-parity implies that the lightest KK mode is stable and can 
be a cold dark matter candidate. 

To begin our calculation we note that the beta functions can be derived more easily in the superficld formalism, 
where we shall now briefly discuss A/" = 1 supersymmetry in a five-dimensional Minkowski s pac e and its description 
in terms of 4D superfields. Note that more details of this approach can be found in Ref s . |15l - [l9l ] . where, as a matter 
of notation, we shall label our space-time coordinates by (x M , y). 

As an introductory example to this formalism, consider the gauge sector, described by a 5D Af — 1 vector super- 
multiplet which consists (on-shell) of a 5D vector field A M , a real scalar S and two gauginos, A and A'. The action 
for which can be given by: 
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S - i\T M D M \ - iX T M D M X' + (A + A )[S, A + A'] 



(2) 
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with Dm = 8m + iAm and T M = (7^, 17 s ). F MN — —^[D M , D N ] and k normalises the trace over the generators of 
the gauge groups. One can also write these fields in terms of a Af = 2, 4D vector supermultiplet, ft — (V , x), where 
V is a A/" = 1 vector supermultiplet containing and A, and x is a chiral Af = I supermultiplet containing A' and 
S' = S + iA 5 . This form of writing the fields follows from the decomposition of the 5D supercharge (which is a Dirac 
spinor) into two Majorana-type supercharges, which constitute a Af = 2 superalgebra in 4D. Both V and x (and their 
component fields) are in the adjoint representation of the gauge group Q. Using the supermultiplets one can write 
the original 5D Af = 1 supersymmetric action Eq.® in terms of Af = 1 4D superfields and the covariant derivative 
in the y direction (l8j : 
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(W a W a S(0 ) + h.c) + (e 



-2gV 



2gV^2 
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with W a 
and V„ = 



= — \D e 2gV D a e 2aV . D a is the covariant derivative in the 4D Af = 1 superspace (see Refs. [20L [2l|) 
d y + X- This action can be expanded and quantised to find the Feynman rules to a given order in the 
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FIG. 1: The one-loop diagrams related to the wave-function renormalisation of the matter superfields, in which diagrams a)-e) 
refer to the case where all the matter Gelds are in the bulk, and the excited KK states are labeled by the number without the 
bracket; whereas diagrams a), c) and d) are related to the brane localised matter fields case, in which the KK states are labeled 
by the number inside the bracket fldl] - 



gauge coupling <?. The details of this procedure can be found in Ref. [lq] . Due to the non-renormalisation theorem 
[221 ] . the beta functions for the couplings of the operators in the superpotential are governed by the wave function 
renormalisation constants Zij = 1 + SZy. The Feynman diagrams related to the wave-function renormalisation are 
given in Fig IT] 

At this point we can note that Higgs superfields and gauge superfields will always propagate into the fifth dimension 
in our model. However, different possibilities of localisation for the matter superfields can be studied by taking the 
two limiting cases of superfields containing the SM fermions in the bulk or all superfields containing SM fermions 
restricted to the brane, respectively. For the case where all fields can propagate in the bulk, the action for the matter 
fields would be flal: 



(4) 



Again, this action can be expanded and quantised, as detailed in Ref. [15|. Similarly, we can do the same for the case 
where all superfields containing SM fermions are restricted to the brane. In which case the part of the action involving 
only gauge and Higgs fields is not modified, whereas the action for the superfields containing the SM fermions becomes 



d s zdyS(y) 2g$ i V$ i } 



(5) 



Before proceeding further we shall recall that the Yukawa couplings in the bulk are forbidden by the 5D M = 1 
supersymmetry. However, they can be introduced on the branes, which are 4D subspaces with reduced supersymmetry 



supei 

m 



■S = / • { ^Ay** { * i <M(0)-h 



(6) 



Returning to FigsfTJ^-e) these lead to the usual Minimal Supersymmetric Standard Model (MSSM) if all the extra- 
dimensional effects were excluded, that is, where the zero modes in FigQJi) and FigUJ:) give the beta functions of the 
usual 4D MSSM. FigfTp) is a new extra-dimensional effect related to the coupling to the chiral superfield originating 
from the fifth component of the gauge fields, and where the higher KK modes in FigsQJi), c), d) and e) are also 
contributing additional extra-dimensional effects. Note that the mass dimension for the Yukawa couplings in 5D is 
— 5, therefore, the power law dependence of their evolution, which results from these radiative corrections, becomes 
significantly important. Also, it should be noted at this point that though the matter fields do not normally receive 
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a wavefunction renormalization in M = 2 supersymmetry, that is the KK conserved diagrams would not contribute 
to the Yukawa evolution, J\f = 2 supersymmetry has been explicitly broken by the compactification [TEj . 

As such, following the procedures outlined in Ref. [l5[ the RGEs for the Yukawa couplings in the 5D MSSM, for 
all three generations propagating in the bulk, is expressed as: 

16^ 2 ^ = Y d (3Tr(YjY d ) + Tr(Y^Y e ) + 3Y^Y d + Y^Y u )irS{tf - Y d (^Lgj + 3g\ + y 5 2 ) S(t) , 

16tt 2 ^ = Yu (3Tr(Y:Y u )+3Y^Y u + YjY d )nS(t) 2 -Y u (^gl + 3gl + fg^S(t), 

16tt 2 ^ = Y e (3Tr(YjY d ) + Tr(Y^Y e ) + 3Y c ?Y e )nS(t) 2 Y e (^gj + 3g^j S(t) . (7) 

That is, when the energy scale /i > 1/R or when the energy scale parameter t = \n(/i/Mz) > ln( M * R ) (where we 
have set Mz as the renormalisation point, and use S(t) = e t MzR) the beta functions become as above, replacing the 
usual 4D MSSM RGEs. Note that, since we compactify with the symmetry y — > —y [S\/Zi symmetry), the integral 
of the extra-dimension is from — >■ irR, not from — > 2ttR. As such the X$ parameter is in our case only half of the 
volume of the 5-dimensional unit sphere (see Ref. 3), i.e. X\ — 1, and this leads to a factor of two difference (due to 
the compactification) with respect to other Refs. [2|. Il5j. 

Furthermore, when < t < ln( ^ R ) (that is, Mz < fi < 1/R) the Yukawa evolution equations at low energy are 
dictated by the usual MSSM, where the above equations reduce to the 4D MSSM equations: 

16^ 2 ^ = Y d (3Tr(YjY d )+Tr(Y:Y e ) + 3YjY d + Y^Y u )-Y d (^gf + 3gi + ^g^ , 

16tt 2 ^ = Y u {3Tr{Y^Y u ) + 3Y^Y U + Y J Yd ) Y u (^g{ + 3g\ + ^ , 

16tt 2 ^ = Y e (3Tr(YjY d ) + Tr(Y e %) + 3Y^Y e ) - Y e (^gj + 3g^j . (8) 

Note also that from the KK number conserved one- loop diagram in FigOJi), at each KK level the new excited states 
exactly mirror the zero mode ground state, where their contributions to the anomalous dimensions are exactly the 
same as those in the usual 4-dimensional MSSM. As a result, when the energy fi — 1/R, or S(t) — 1, the part in the 
beta function associated with the gauge fields in Eqs.© reduce to the normal 4D MSSM formalism in Eqs.(|B])- 

For the second case we shall consider, that of the matter fields localised on the brane, we have (in the Feynman 
diagrams, only FigOJ,), c) and d) will contribute to the wave function of the matter fields, while the non-zero integer 
numbers inside the bracket labeling the KK states of the gauge and Higgs superfields respectively) 

167r2 ^T = Y ^ Tr ( Y l Y ^ + Tr ( Y e Y e) + (GYjY d + 2Y^Y u )S{t)) Y d {^g\ + \g\ + |.g 3 2 ) S(t) , 
16^ 2 ^ = Y u (3Tr(Y^Y u ) + (QY^Y U + 2YjY d )S{t)) - Y u (^g\ + \g\ + S(t) , 

16tt 2 ^ = r e (3Tr(Yjr d )+Tr(r;r e ) + ^^^^ 

As such, if we write the evolution of the gauge couplings in 4-dimensions as: 

16tt 2 ^* = b l9 ? , (10) 

where in the 4D MSSM the parameters hi read (&i, 62,^3) = 3) [7], using an SU(5) normalisation. Then if 

we consider our 5D theory as effective up to a scale A, then between the scale i? -1 (the compactification scale of the 
single flat extra-dimension), where the first KK states are excited and the cut-off scale A, for the gauge couplings, 
there are finite quantum corrections from the AR number of KK states. As a result, once the KK states are excited, 
these couplings exhibit power law dependencies on A. This can be illustrated if AR > 1, to a very good accuracy, 
with the generic 4D beta function of the gauge couplings with the power law evolution behaviour [4j, [5j 

(3 iD ^p iD + (S(»)-1)P, (11) 



(9) 
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FIG. 2: (Colour online) Gauge couplings (g\ (red), g2 (blue), g^ (green) with: in the left panel, all matter fields in the bulk; 
and the right panel for all matter fields on the brane; for three different values of the compactification scales (2 TeV (solid 
line), 8 TeV (dot-dashed line), 15 TeV (dashed line)) as a function of the scale parameter t. 



where /? is a generic contribution from a single KK level, and where its coefficient is not a constant but instead 
S(p) = /iR, with fj, max = A, reflecting the power law running behaviour. Therefore, in terms of the scale parameter 
t, the evolution of the gauge couplings can be written as: 

I6vr 2 ^ = [h + (S(t) - l)bi] gi 3 . (12) 

Next we consider the beta functions of the gauge couplings in the 5D MSSM. In fact, after compactification of the 
5D MSSM, we have two 4D Af = 1 chiral supermultiplets, $ and $ c , where the zero modes of $ give us the normal 
matter fields and Higgs fields as well as their super partners, while $ c is a new supermultiplet. In the component 
field formalism, at each KK level, aside from the quantum corrections that mirror those of the 4D MSSM, the only 
new one-loop contributions to the Feynman diagrams are from the wave function renormalisation of (which 
contribute via the coupling of A^ to the complex scalar field and its super-partner in the superfield x, an d the coupling 
of A^ with the new fermion field and its super-partner in the superfield $ c associated with the two doublets of the 
Higgs fields and the matter fields respectively in the bulk). This then gives rise to the master beta functions of the 
gauge couplings in the 5D MSSM as follows: 

(61,62,63) = 0,-2,-6^+477, (13) 

where 77 represents the number of generations of fcrmions which propagate in the bulk. So in the two cases we shall 
consider, that of all fields propagating in the bulk, i.e. rj — 3, we have |4|: 

bi = (™ ,10,6) . (14) 

Similarly, for all our matter fields localised to the 3-brane (that is, r\ — 0), we have: 

bi = Q,-2,-V) . (15) 

In FigJ5]we have plotted the evolutions of the brane localised and bulk field cases for several choices of compact- 
ification scales for the extra-dimension (R). From these plots, and the discussion given in Ref.^], we find that for 
the three gauge coupling constants to approach a small region at some value of t requires an extremely large value 
of 1/R, which is of no phenomenological interest at present. For the case of our fields being brane localised, the 
extra-dimensions naturally lead to gauge coupling unification at an intermediate mass scale for the compactification 
radii considered here. 
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III. SCALING OF THE CKM MATRIX 



The RGEs are an important tool for investigating the properties of the quark masses and the CKM matrix at 
different energy scales, this therefore strongly constrains the possible symmetries or textures at the grand unification 
scale. Having derived a set of equations for the evolutions of the Yukawa and gauge couplings in the previous section, 
we now develop the evolutions of the CKM matrix in the 5D MSSM for our two cases. Consider first the case of all 
the matter fields propagating in the bulk. From the Yukawa evolution equations Eqs.([7]) we can convert them to the 
following form: 

,dY d 



16tt 2 ^ = Y d {T d irS 2 -G d + (3Y}Y d + YjY u )TrS 2 } , 



16^ 2 ^ = Y u {t u ttS 2 -G u + (3YJY u + Y}Y d )irS 2 } , 
dY 

16^ 2 ^ = Y e {T e nS 2 ~G e + {3Y^Y e )^S 2 } , (16) 

where T d = 3Tr{Y^Y d ) + Tr(Y^Y e ), G d = (^g 2 + 3g 2 + fg 2 ) S(t), T u = 3Tr(Y^Y u ),G u = (±§ g\ + 3g\ + fg 2 ) S(t), 

T e = 3Tr(YjY d ) + Tr(Y^Y e ) and G e = (|g 2 + 3g\) S(t). Further, the evolution of the square of the Yukawa coupling 
matrices becomes: 

167T 2 ^J^ = 2(T d nS 2 - G d )(YjY d ) + 6(YjY d ) 2 nS 2 + [(Y^Y u )(Y^Y d ) + {Y^ d Y d ){Y^Y u )]^S 2 , 
16 ^^T^ = 2{TunS2 ~ G ^ Y ^ + 6 ( Y u Y ^ s2 + [{y^u){Y ] d Y d ) + {Y}Y d ){Y^Y u )]vS 2 , 

^ ^dt^ = 2 ( T ^ 2 ~ G °)( Y 2 Y ') + HY e ^Y e ) 2 nS 2 . (17) 
The square of the quark Yukawa coupling matrices can be diagonalized by using two unitary matrices U and V, 

diag(/ 2 ,/ 2 ,/ 2 ) - UYfcrf, 

diag (h 2 ,h 2 s ,h 2 ) = VYjY d V^, (18) 

in which f 2 , f 2 , f 2 and h 2 d , h 2 , hi are the eigenvalues of Y^Y U and Y d Y d respectively. It follows that the CKM 
matrix appears as a result of the transition from the quark flavor eigenstates to the quark mass eigenstates upon this 
diagonalization of the quark mass matrices: 

Vckm = UV^ . (19) 

By imposing the unitary transformation Ea. (TT8|) . on both sides of the evolution equations of Y^Y U and Y d Y d , and 
taking the diagonal elements, we obtain the following two relations: 

Wn 2 ^- = f 2 [2(T u nS 2 -G u )+6nS 2 f 2 + 2nS 2 Y,h 2 \V tJ \ 2 } 

3 

dh 2 

16tt 2 ^ = h 2 [2(T d nS 2 -G d ) + 6nS 2 h 2 + 2nS 2 Y,f?\^j\ 2 } , (20) 

i 

where Vij are the elements of Vckm- These along with the following equation 

dn 2 

167 r2 ~Jr = y ^ TeWS ~ Ge) + 6nS2y2] ' (21) 

for the lepton sector, constitute a complete set of coupled differential equations for the three families, in which 
y 2 e = diag (y 2 , y 2 , y 2 ) = Y^Y e . 

Considering the variation of the square of the quark Yukawa couplings, we may impose two new unitary matrices to 
make them diagonal. Thus, by applying Eq. fp~9|) . we are led to the variation of the CKM matrix and thus its evolution 
equation when the energy scale is beyond the threshold R . 



. a 2 dVik c2 
167T — - — = nb 
dt 



f 2 + f 2 h 2 + h 2 

72 _ j2 h m V imVj* m V jk + ^ h 2 _ j£fj V i*m V jk V ir 



(22) 
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Likewise, for the case where the matter fields are localised to the brane, Eq.© gives us the following: 

,dY d 



16tt 2 ^ = Y d \T d - G d + (6YjY d + 2Y^Y U )S] 



16tt 2 ^ = Y u {T u -G u + {6Y^Y u + 2YjY d )s} , 
iTY 

16 ^^f = Y e {T e ~G e + (6Y^Y e )S} , 



(23) 



where T d = 3Tr(YjY d ) + Tr(Y}Y e ), G d = (|§ 5 2 + § 5 2 + f g 2 ) S(t), T u = 3TrYjY a , G u = {§g 2 + \g\ + fg 2 ) S(t), 
T e = 3Tr(YjY d ) + Tr(Y^Y e ) and G e = (f§ 

9i + ffli) ^(t). As was done earlier, these can be manipulated to: 



16 

16tt 



16. 



: 2{T U - G U ){Y^Y U ) + \2{Y^Y U ) 2 S + 2[(Y^Y u )(YjY d ) + (Y^Y d )(Y^Y u )]S , 

: 2(r e -G e )(y e ty e ) + i2(y e ty e ) 2 s . 



dt 

2 d(YjY u 
dt 

2 d(YjY e 



dt 



(24) 



Giving us 



16.* = ffWTu-G^ + USfi + iSY^tflVitf] 

3 

dh 2 

16tt 2 ^ = h 2 [2(T d -G d ) + l2Sh 2 + 4Sj2f?\Vv\ 2 ]-. 

i 

y 2 e [2(T e -G e ) + 12Sy 2 e ] , 



le. 2 ^ 



dt 

as well as the evolution equation of the CKM matrix elements 



16.^ 



dt 



= 25 



f 2 4- f 2 h 2 A- h 2 

E J J___Jj_h 2 V. V* V J- \^ k m ? 2 ^r* T/ t/. 



f 2 - f 2 



(25) 



(26) 



in which the ff, h 2 , y\ are the eigenvalues of Y^Y u ,Y^Y d , Y^Y e as defined earlier. 

However, when the energy is below the compactification scale, we have the usual 4D MSSM evolution equations for 
the Yukawa couplings and the CKM matrix 



16.* 
dt 


= f?[2(T u 


-G„)- 


f 6/ 2 - 


j 


? dh 2 
16. 2 — ± 
dt 


= h 2 [2(T d 


-G d )- 




-sE/.i^i 2 ] 

■j 


16.* 
dt 




-G e )4 


- 6y 2 ] , 





and 



16. 



dt 



/? + / 



hl + hl 



El±_JJ_h 2 V- V* V J- "fc "m j2t/-* t/ t/. 

f 2 _ f2 "m ™ K jra 'jk 1" j,2 f,2 *j v jm v jk »im 



f? - f 2 



J, 2 _ U2 ->3 • 3-r 



(27) 



(28) 



where T u = 3Tr(y„tF u ), G„ = if.? 2 + 3<? 2 + fg 2 , T d - 3Tr(rjr d ) + Tr(y e tr e ), G d = ^.g 2 + 3 5 2 + f 3 2 , T e = 

3Tr(y]r d ) + Tr{Y}Y e ) and G e = f g 2 + 3g 2 . 

Before proceeding to numerically solve these systems of equations, we choose the initial input values for the gauge 
couplings, fermion masses and CKM elements at the Mz scale (as in Ref.[5]). The 4D MSSM contains the particle 
spectrum of a two-Higgs doublet model extension of the SM and the corresponding supersymmetric partners. After 
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Runnings of /, for tan/?= 1 



Runnings of /, for tan/3=30 
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FIG. 3: (Colour online) The Yukawa coupling ft for the top quark (in the bulk) as a function of the scale parameter t, for 
(left panel) tan/3 = 1 and (right panel) tan/3 = 30 for different compactification scales: R" 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 



the spontaneous breaking of the electroweak symmetry, five physical Higgs particles are left in the spectrum. The 
two Higgs doublets H u and H d , with opposite hypercharges, are responsible for the generation of the up-type and 
down-type quarks respectively. The vacuum expectation values of the neutral components of the two Higgs fields 

satisfy the relation v u 2 + V4 2 — (^7§) = {^7AGeV) 2 . The fermion mass matrices appear after the spontaneous 
symmetry breaking from the fermion-Higgs Yukawa couplings. As a result, the initial Yukawa couplings are given by 
the ratios of the fermion masses to the appropriate Higgs vacuum expectation values as follows: 

r <^u,c,t 7 ^d,s,b ^e,/_i,r 

Ju,c,t — j ^d,s,6 — 5 Ue,iiT 7 yZu J 

v u v d v d 

where we define tan /3 = v u /vd, which is the ratio of vacuum expectation values of the two Higgs fields H M and H d . 



IV. DISCUSSION OF THE RESULTS 



In the previous section we have derived the full set of one-loop coupled RGEs for the Yukawa and gauge couplings, 
together with the CKM matrix elements for both the universal 5D MSSM and brane localised matter field scenarios. 
From these complete sets of the RGEs we can obtain the renormalisation group flow of all observables related to up- 
and down-quark masses and the quark flavor mixings. For our numerical analysis we assume the fundamental scale 
is not far from the range of LHC scale, and set the compactification scale to be = 2 TeV, 8 TeV, and 15 TeV 
respectively. 

Actually, below the supersymmetric breaking scale the Yukawa and gauge couplings run in the usual logarithmic 
fashion, giving a rather slow change for their values. Therefore, for supersymmetric breaking theories around TeV 
scales, for simplicity, we take the supersymmetric breaking scale Msusy = Mz in the present numerical study, and 
run the RGEs from Mz up to the high energy scales for our three different compactification scales. In Fig(5]we find 
the evolution of the gauge couplings have different properties for these two 5D MSSM scenarios. 

Additionally, as illustrated in FigJ31 for the case of universal 5D MSSM, once the first KK threshold is crossed at 
fx = R , the power law running of the various beta functions causes the Yukawa coupling to rapidly increase following 
the rapid increase in the gauge coupling constants in the left panel of FigJ5] From Eq. (|2T)|) we can find the quadratic 
term of S(t) providing a positive contribution to the Yukawa beta functions, which is in contrast to beta functions 
of the gauge couplings (which include terms only linear in S(t)). Therefore, from Eq.(|20j). the positive contribution 
from S(t) terms will dominate the negative contributions from the gauge couplings, and cause the Yukawa couplings 
to increase rapidly. This behaviour can be observed for both small and large tan j3 cases. However, as illustrated 
in the first graph of Fig|3l for small tan/3, the Yukawa coupling has a large initial value, therefore it blows up at a 
relatively low energy as compared with the case for large tan/3. As a result, as one evolves upward in the scale, the 
top Yukawa coupling is rising with a fast rate and is pushed up against the Landau pole where it becomes divergent 
and blows up. In the vicinity of this singular point the perturbative calculation becomes invalid, and the higher order 
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FIG. 4: (Colour online) The Yukawa coupling ft for the top quark (on the brane) as a function of the scale parameter t, for 
(left panel) tan/3 = 1 and (right panel) tan/3 = 30 for different compactihcation scales: R" 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 



corrections become significant. The Landau pole also indicates that there is an upper limit on the value of the gauge 
couplings where new physics must emerge before the Yukawa couplings diverge. 

In the brane localised matter field scenario, the beta function has only linear terms in S(t), which is comparable 
with the S(t) term in the beta function for the gauge couplings. As depicted in FigJH for a small value of tan/3, 
we have a large initial value of ft and the gauge coupling contribution to the Yukawa beta function is sub-dominant 
only. Therefore, as discussed previously, the Yukawa coupling f t increases rapidly as one crosses the KK threshold 
at /i = i? -1 , resulting in a rapid approach of the singularity before the unification scale is reached. However, for an 
intermediate value of tan /3, we have a relative smaller initial condition for the top Yukawa coupling and the Yukawa 
terms in the beta function become less important. The contributions from the gauge couplings may then become 
significant, which leads to a net negative contribution to the beta functions. Therefore, the curvature of the trajectory 
of the top Yukawa evolution might change direction, and the Yukawa evolution will decrease instead of increasing. 
This behaviour would become more obvious for a large value of tan /3. As observed in FigJU for tan /3 — 30, we indeed 
observe the decreasing behaviour of the top Yukawa couplings. This behaviour provides a very clear phenomenological 
signature, especially for scenarios with a larger tan /3 and that are valid up to the unification scale where the gauge 
couplings converge. 

We next turn our attention to the quark flavor mixing matrix, especially the complex phase of the CKM matrix 
which characterises CP-violating phenomena. This phenomena has been unambiguously verified in a number of K — K 
and B — B systems. Because of the arbitrariness in choice of phases of the quark fields, the phases of individual matrix 
elements of the Vqkm are not themselves directly observable. Among these we therefore use the absolute values of 
the matrix element \Vij\ as the independent set of rephasing invariant variables. Of the nine elements of the CKM 
matrix, only four of them are independent, which is consistent with the four independent variables of the standard 
parametrisation of the CKM matrix. For definiteness we choose the | V u b | , | V c b \ , | V us | and the Jarlskog rephasing 
invariant parameter J = lmV u dV cs V* s V* d as the four independent parameters of Vckm- In FigslSJH] we plot the 
energy dependence of these four variables from the weak scale all the way up to the high energy scales for different 
values of compactification radii R^ 1 for the universal 5D MSSM case, and in FigsJ7][5]for the brane localized matter 
fields case. In these sets of pictures we consider two indicative choices of tan (3, that of tan (3 — 1 and tan /3 = 30. 

The running of the CKM matrix is governed by the terms related to the Yukawa couplings. These Yukawa couplings 
are usually very small, except for the top Yukawa coupling (which could give a sizeable contribution). The CKM 
matrix elements V us ~ 6*i 2 , V u b — 9i3e~ lS , V c b — #23, can be used to observe the Cabibbo angle, #13 and #23- For the 
universal 5D MSSM scenario we plot their evolution in FigslSJH and find that they decrease with the energy scale in 
a similar manner regardless of whether tan/3 is small or large. However, for a large initial value of ft (small tan/3), 
the mixing angles have a more rapid evolution and end in the regime where the top Yukawa diverges and develops 
a singularity. Similarly, the Jarlskog parameter also decreases quite rapidly once the initial KK threshold is passed. 
However, when tan /3 is large, we have a relatively longer distance between the initial and terminating energy track, 
the evolution of J can be driven towards zero or even further. Besides, as can be seen explicitly in Ref. Q, the beta 
functions of the evolution equations of the CKM elements are up to the third order in the CKM elements, which are 
comparably smaller than that of Jarlskog parameter's quadratic dependence on the CKM elements. This fact then 
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FIG. 5: (Colour online) The CMK matrix elements \V US \ (top left), \V u b\ (top right), \V c b\ (bottom left) and the Jarlskog 
parameter J (bottom right) as functions of the scale parameter t for tan /3 — 1. All matter fields are in the bulk for a variety 
of compactification scales: R' 1 = 2 TeV (red, dotted line), 8 TeV (blue, dot-dashed line), and 15 TeV (green, dashed line). 



leads to the relatively large variation of J with the increase of energy. Quantitatively we observe from these plots the 
following; the decrease in the values of V us is not sizeable, whilst for V u b and V c b their variations change by more than 
50%. Furthermore, for tan/3 = 30, the Jarlskog parameter drops almost to zero, which sets the effect of the SM CP 
violation to being very small. Note, however, that in a supersymmetric theory other sources of CP violation beyond 
the SM ones are typically present, therefore only a complete and detailed study of a specific model would allow us to 
establish the strength of the CP violating effects. 

For the case of the matter fields being constrained to the brane, in Figsj7][8]we observe that the evolutions of these 
mixing angles and CP violation parameter are decreasing irrespective of whether the top Yukawa coupling grows or 
not. For small tan/3 we see similar evolution behaviours for these parameters as in the bulk case. The decreases in 
the values of V us , V u b, V c b and J are much steeper, due to rapid growth of the top Yukawa coupling near the singular 
point. However, as tan/3 becomes larger, e.g. tan/3 = 30, the top Yukawa coupling evolves downward instead of 
upward. The decreases in these CKM parameters then becomes much milder towards the unification scale; though 
the reduction to effectively zero in the Jarlskog parameter persists. As a result, for the brane localised matter field 
scenario, it is more desirable to have a large tan (3 for theories that are valid up to the gauge coupling unification 
scale. 



V. CONCLUSIONS 



In summary, for the two 5D MSSM scenarios with matter fields in the bulk or on the brane, we have performed the 
numerical analysis of the evolution of the various parameters of the CKM matrix, and both cases give us a scenario 
with small or no quark flavor mixings at high energies, especially for the mixings with the heavy generation. The 
evolution equations which relate various observables at different energies, and also allow the study of their asymptotic 
behaviours, are particularly important in view of testing the evolution of the Yukawa couplings. In the universal 5D 
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FIG. 6: (Colour online) The CMK matrix elements \V US \ (top left), \V u b\ (top right), \V c t\ (bottom left) and the Jarlskog 
parameter J (bottom right) as functions of the scale parameter t for tan f) — 30. All matter fields are in the bulk for a variety 
of compactification scales: R~ x = 2 TeV (red, dotted line), 8 TeV (blue, dot-dashed line), and 15 TeV (green, dashed line). 



MSSM model, the evolution of these CKM parameters have a rapid variation prior to reaching a cut-off scale where 
the top Yukawa coupling develops a singularity point and the model breaks down. For the brane localised matter 
fields model, we can only observe similar behaviours for small values of tan/?, while for large tan/3, the initial top 
Yukawa coupling becomes smaller, the gauge couplings then play a dominant role during the evolution of the Yukawa 
couplings, which cause the Yukawa couplings to decrease instead of increasing. As such the variations of these CKM 
parameters have a relatively milder behaviour, and the theory is valid up the gauge coupling unification scale. 
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FIG. 8: (Colour online) The CMK matrix elements \V US \ (top left), \V u b\ (top right), \V c b\ (bottom left) and the Jarlskog 
parameter J (bottom right) as functions of the scale parameter t for tan ft = 30. All matter fields are on the brane for a variety 
of compacti&cation scales: R -1 = 2 TeV (red, dotted line), 8 TeV (blue, dot-dashed line), and 15 TeV (green, dashed line). 



